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Abstract  To  achieve  a  tight  integration  of  design  and  anal¬ 
ysis,  conformal  solid  T-spline  construction  with  the  input 
boundary  spline  representation  preserved  is  desirable.  How¬ 
ever,  to  the  best  of  our  knowledge,  this  is  still  an  open  prob¬ 
lem.  In  this  paper,  we  provide  its  first  solution.  The  input 
boundary  T-spline  surface  has  genus-zero  topology  and  only 
contains  eight  extraordinary  nodes,  with  an  isoparametric 
line  connecting  each  pair.  One  cube  is  adopted  as  the  para¬ 
metric  domain  for  the  solid  T-spline.  Starting  from  the  cube 
with  all  the  nodes  on  the  input  surface  as  T-junctions,  we 
adaptively  subdivide  the  domain  based  on  the  octree  struc¬ 
ture  until  each  face  or  edge  contains  at  most  one  face  T- 
junction  or  one  edge  T-junction.  Next,  we  insert  two  bound¬ 
ary  layers  between  the  input  T-spline  surface  and  the  bound¬ 
ary  of  the  subdivision  result.  Finally,  knot  intervals  are  cal¬ 
culated  from  the  T-mesh  and  the  solid  T-spline  is  constructed. 
The  obtained  T-spline  is  conformal  to  the  input  T-spline  sur¬ 
face  with  exactly  the  same  boundary  representation  and  con¬ 
tinuity.  For  the  interior  region,  the  continuity  is  C2  every¬ 
where  except  for  the  local  region  surrounding  irregular  nodes. 
Several  examples  are  presented  to  demonstrate  the  perfor¬ 
mance  of  the  algorithm. 
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1  Introduction 

Spline  boundary  representation  is  the  most  popular  technol¬ 
ogy  for  solid  modeling  in  computer  aided  design  (CAD). 
To  integrate  design  and  analysis,  the  root  idea  of  isogeo¬ 
metric  analysis  [6,  2],  one  challenge  is  to  automatically  cre¬ 
ate  a  conformal  solid  NURBS/T- spline  model  with  the  given 
spline  boundary  representation  preserved  exactly.  To  the  best 
of  our  knowledge,  this  paper  presents  the  first  attempt  to 
construct  conformal  solid  T- splines  from  boundary  spline 
representations. 

The  advent  of  T-splines  [10,  9]  provides  not  only  more 
flexibility  for  modeling,  but  also  an  ideal  discretization  for 
isogeometric  analysis  [8].  Several  methods  have  been  devel¬ 
oped  for  solid  spline  modeling.  Hua  et  al.  [5]  developed  a 
trivariate  simplex  spline  model  which  is  defined  over  a  tetra¬ 
hedral  decomposition  of  any  3D  domain  with  complicated 
geometry  and  arbitrary  topology.  In  [13],  a  skeleton-based 
solid  NURBS  construction  method  for  patient- specific  vas¬ 
cular  geometric  models  was  presented.  In  [1],  a  swept  vol¬ 
ume  parameterization  was  built  for  isogeometric  analysis.  A 
framework  was  developed  in  [7]  to  model  a  single  trivariate 
B- spline  from  input  boundary  triangle  meshes  with  genus- 
zero  topology.  Escobar  et  al.  [4]  proposed  a  solid  T-spline 
modeling  algorithm  based  on  optimization  from  the  surface 
triangulation.  In  our  earlier  work,  we  developed  rational  T- 
splines  whose  basis  functions  satisfy  partition  of  unity  by 
definition  [12].  In  addition,  we  developed  solid  T-spline  con¬ 
struction  methods  for  both  genus-zero  geometry  [14]  and 
high  genus  geometry  [11]  from  boundary  triangle  meshes, 
using  a  unit  cube  or  a  polycube  as  the  parametric  domain. 
These  developments  provide  a  foundation  upon  which  the 
conformal  trivariate  spline  model  can  be  built  from  spline 
boundary  surface  representations. 

A  general  methodology  for  constructing  a  conformal  solid 
T-spline  from  boundary  T- spline/NURBS  representations  is 
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Fig.  1  Flow  chart  of  conformal  solid  T-spline  construction  from  given  boundary  NURBS/T- spline  surface  representations  with  genus-zero  topology 
and  eight  extraordinary  nodes,  (a)  The  input  boundary;  (b)  octree  subdivision  &  projection;  (c)  first  boundary  layer  construction;  (d)  second 
boundary  layer  construction;  (e)  handling  irregular  nodes;  and  (f)  the  constructed  solid  T-spline  with  T-mesh. 


very  complex  due  to  complicated  geometry  and  topology, 
various  T-mesh  or  NURBS  patch  connectivity,  and  various 
parametrizations.  The  conformal  property  is  hard  to  main¬ 
tain  during  the  process  of  converting  a  bivariate  surface  model 
to  a  trivariate  solid  model.  Here,  we  focus  on  geometry  with 
genus-zero  topology.  In  addition,  the  input  T-spline  surface 
only  contains  eight  extraordinary  points. 

In  this  paper,  we  propose  a  novel  method  to  construct 
conformal  solid  rational  T- splines  for  genus-zero  geometry, 
with  the  input  boundary  T-spline  representation  preserved 
exactly.  We  use  one  cube  as  the  parametric  domain  for  the 
solid  T-spline.  Also,  an  adaptive  subdivision  is  applied  to 
the  cube  and  the  boundary  nodes  obtained  in  the  subdivi¬ 
sion  are  projected  to  the  input  T-spline  surface  via  the  sur¬ 
face  parametrization.  To  conform  to  the  input  boundary,  two 
boundary  layers  are  generated  between  the  input  boundary 
surface  and  the  boundary  of  the  subdivision  result.  Tem¬ 
plates  for  irregular  nodes  are  then  applied  to  get  a  gap-free 
T-mesh.  Finally,  a  solid  rational  T-spline  is  built  from  the 
T-mesh  and  solid  Bezier  elements  are  extracted. 


The  remainder  of  this  paper  is  organized  as  follows.  Sec¬ 
tion  2  presents  an  overview  of  the  method.  The  octree  subdi¬ 
vision  algorithm  is  described  in  Section  3  and  the  boundary 
layer  method  is  explained  in  Section  4.  Section  5  discusses 
the  method  to  deal  with  irregular  nodes  and  solid  T-spline 
construction.  Section  6  presents  some  examples  and  Section 
7  draws  conclusions. 


2  Problem  Description  and  Algorithm  Overview 

Figure  1  shows  the  pipeline  of  conformal  solid  T-spline  con¬ 
struction.  The  input  can  be  a  regular  T-spline  surface, 


£  vv,C, 5,(^1!) 
S&Ti)  =  ^r - , 

|) 

i= 0 


(5.ti)  en, 


(i) 


Conformal  Solid  T-spline  Construction  from  Boundary  T-spline  Representations 


3 


(d)  (e)  (f) 

Fig.  2  The  subdivision  process  for  one  2D  example,  (a)  The  input  boundary  curve  with  its  control  polygon  in  the  physical  space;  (b)  subdivision 
result  after  one  refinement;  and  (c)  final  subdivision  result,  (d-f)  shows  the  corresponding  result  in  the  parametric  space  for  (a-c),  respectively. 


or  a  rational  T-spline  surface, 

- ,  Men,  (2) 

i= 0 


where  Wj  is  the  weight  for  the  control  point  Q,  |)  = 
^W(Tl),AfandA^  are  B -spline  basis  functions  defined 
by  two  local  knot  vectors  at  node  Q,  7^  (^,r|)  is  the  rational 


B -spline  basis  function,  R;(^,r|)  = 


^=0ivf(^(Tl) 


,  and 


Q  is  the  domain  of  the  T-spline  surface.  For  our  conformal 
solid  T-spline  construction,  the  input  surface  should  satisfy 
the  following  requirements: 


-  Its  topology  has  to  be  genus-zero,  and  we  use  one  cube 
C  as  the  parametric  domain  of  the  solid  T-spline. 

-  The  input  surface  must  contain  eight  extraordinary  nodes 
Q  (i  =  0, . . . ,  7)  and  there  is  an  isoparametric  line  be¬ 
tween  each  pair  of  extraordinary  nodes  QCy,  Extraordi¬ 
nary  nodes  in  2D  are  irregular  ones  with  a  valence  other 
then  four.  Each  extraordianry  node  here  is  mapped  to 
one  corner  of  the  cube  C- 

-  By  splitting  the  input  surface  using  the  isoparametric 
lines  QCj,  each  face  of  the  cube  C  can  serve  as  the  para¬ 
metric  domain  for  its  corresponding  patch. 


In  summary,  the  boundary  of  the  cube  can  serve  as  the  para¬ 
metric  domain  of  the  input  boundary  surface. 

As  shown  in  Figure  1,  given  the  input  NURBS/T- spline 
surface  shown  in  (a),  we  first  use  octree  subdivision  to  ob¬ 
tain  an  initial  T-mesh  shown  in  (b).  Then  two  boundary  lay¬ 


ers  shown  in  (c)  &  (d)  are  built  in  order  to  preserve  the  in¬ 
put  boundary  representation.  Then,  we  deal  with  irregular 
nodes  and  also  insert  zero-knot-interval  edges  to  get  a  valid 
T-mesh.  Finally,  we  calculate  the  knot  vectors  for  each  node 
and  construct  trivariate  solid  T- splines. 


3  Octree  Subdivision  and  Projection 

We  use  a  cube  C  as  the  parametric  domain  of  the  constructed 
solid  T-spline.  The  eight  extraordinary  nodes  correspond  to 
the  eight  comers  of  the  cube.  The  twelve  isoparametric  curves 
connecting  these  extraordinary  nodes  are  associated  with  the 
edges  of  the  cube.  Based  on  these  curves,  we  divide  the  input 
surface  into  six  patches,  and  each  one  is  associated  with  one 
face  of  C-  For  the  input  NURBS/T- spline  surface,  the  para¬ 
metric  domain  of  each  NURBS/T- spline  patch  corresponds 
to  one  face  of  the  cube  C-  Then  we  build  an  initial  T-mesh 
by  applying  the  adaptive  octree  subdivision. 

We  first  assign  all  the  nodes  in  the  input  T-spline  sur¬ 
face  except  the  extraordinary  nodes  as  T-junctions  for  the 
unit  cube  C.  Starting  from  C  with  T-junctions,  each  ele¬ 
ment  is  refined  recursively  into  eight  until  each  boundary 
element  contains  at  most  one  face  T-junction  on  a  face,  and 
at  most  one  edge  T-junction  on  an  edge.  Two  nodes  with 
the  same  parametric  value  are  treated  as  one.  For  solid  T- 
splines,  one  edge  T-junction  is  a  T-junction  lying  on  one 
edge,  and  one  face  T-junction  is  a  T-junction  lying  on  one 
face.  For  each  refinement,  instead  of  always  using  the  cen¬ 
tral  parameter  value,  one  parameter  value  of  the  T-junctions 
in  this  element  is  chosen  to  refine  the  element.  If  there  are 
several  T-junction  parametric  values  in  one  direction,  the 
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one  closest  to  the  middle  is  chosen.  In  this  way,  we  try  to 
minimize  the  number  of  T-junctions.  The  strongly -balanced 
octree  is  adopted  here.  During  subdivision,  the  newly  gen¬ 
erated  boundary  nodes  are  projected  to  their  corresponding 
positions  on  the  input  T- spline  surface.  The  interior  nodes 
are  relocated  to  the  mass  center  of  all  neighboring  elements. 
Then  an  initial  T-mesh  is  generated. 

Figure  2  demonstrates  the  subdivision  process  using  a 
2D  example,  (a)  is  the  input  B -spline  boundary  curve.  The 
dark  blue  polygon  shows  the  boundary  control  mesh  and  the 
light  blue  curves  are  the  B- spline  boundary,  which  has  four 
sharp  comers  rendered  in  red.  The  values  labelled  around  the 
edges  show  the  knot  intervals,  (d)  shows  the  parametric  do¬ 
main  for  the  solid  T-spline,  and  we  assign  nodes  in  the  input 
curve  as  T-junctions  for  the  parametric  domain,  (b)  and  (e) 
are  the  first  refinement  results  in  the  physical  and  parametric 
space,  respectively,  in  which  the  parametric  values  of  Ce  and 
Ci 9  are  used  for  £  and  r|  directions  since  they  are  closest  to 
the  parametric  center  for  this  element.  Then  we  find  that  the 
right  upper  element  contains  two  T-junctions  in  the  £  direc¬ 
tion  (Ci 4  and  C15),  hence  we  subdivide  that  element  again. 
This  process  continues  until  we  obtain  the  result  shown  in 
(c)  and  (f),  for  which  each  boundary  element  contains  only 
one  T-junction  in  each  direction  except  for  the  nodes  which 
have  the  same  parametric  position  with  the  comer  nodes  (Co, 
Ci,  C2  and  C3). 

Figure  3  shows  the  subdivision  result  for  a  Mouse  model, 
(a)  and  (b)  show  the  input  T-spline  surface  in  the  physical 
space  and  the  parametric  space,  respectively,  (c)  and  (d)  are 
the  subdivision  results. 


Fig.  3  The  subdivision  result  for  the  Mouse  model,  (a)  The  input  T- 
spline  surface;  (b)  the  parametric  mapping  result;  (c)  the  subdivision 
result  in  the  physical  space;  and  (d)  the  subdivision  result  in  the  para¬ 
metric  space. 


4  Boundary  Layer  Construction 

We  insert  two  boundary  layers  between  the  subdivision  re¬ 
sult  and  the  input  boundary  surface  in  order  to  make  the 
solid  T-spline  conformal  to  the  input  surface  and  also  im¬ 
prove  the  quality  of  the  T-mesh.  Figure  4  illustrates  the  bound¬ 
ary  layer  constmction  using  the  previous  2D  example,  in 
which  (a)  is  the  subdivision  result,  (b)  is  the  boundary  layer 
result  for  the  subdivision  boundary,  and  (c)  is  the  boundary 
layer  result  for  the  input  boundary  surface.  First,  we  dupli¬ 
cate  each  boundary  node  lying  on  the  element  corner,  insert 
one  element  for  each  boundary  face,  and  obtain  one  bound¬ 
ary  layer  (rendered  in  blue  in  (b)).  We  assign  the  bound¬ 
ary  T-junctions  to  the  newly  generated  boundary  faces.  As 
shown  in  (b),  the  boundary  T-junctions,  such  as  C5,  C7,  C10, 
Ci 5  and  Ci6,  are  assigned  to  the  new  boundary  after  bound¬ 
ary  layer  construction.  Second,  we  duplicate  the  boundary 
nodes  which  have  the  same  parametric  position  with  one  of 
the  input  control  points,  such  as  C5,  C6  and  C7.  The  physical 
position  is  assigned  to  be  the  same  as  the  input.  The  sec¬ 
ond  pillowed  layer  is  rendered  in  yellow  in  (c).  Next,  we  use 
smoothing  and  optimization  techniques  to  relocate  the  phys¬ 
ical  positions  of  the  interior  nodes  to  improve  the  T-mesh 
quality. 

Figure  5  shows  the  result  after  boundary  layer  construc¬ 
tion  for  the  Mouse  model.  In  (b),  the  pink  region  is  the  sub¬ 
division  part,  the  blue  layer  is  the  boundary  layer  for  the 
subdivision  result,  and  the  yellow  layer  is  the  second  bound¬ 
ary  layer.  These  two  layers  serve  as  a  transition  region  be¬ 
tween  the  input  boundary  spline  surface  and  the  subdivision 
result.  From  this  figure,  we  can  observe  that  the  extraordi¬ 
nary  nodes  lie  in  the  interior  and  have  two  layers  of  ele¬ 
ments  away  from  the  boundary  and  we  generated  some  T- 
j unctions  between  the  boundary  layers.  We  can  guarantee 
that  after  this  boundary  layer  construction  step,  the  boundary 
of  the  solid  T-mesh  is  exactly  the  same  as  the  input.  In  other 
words,  the  constructed  solid  T-spline  is  conformal  to  the  in¬ 
put  NURBS/T- spline  surface.  We  formulize  these  statements 
in  a  lemma. 

Lemma  1  Due  to  the  boundary  layer  construction,  the  ob¬ 
tained  solid  T-spline  is  conformal  to  the  input  NURBS/T- 
spline  surface.  In  other  words,  the  input  NURBS/T- spline 
surface  is  preserved  exactly. 

Proof:  Due  to  the  boundary  layer  construction,  the  second 
boundary  layer  can  be  treated  as  one  hexahedral  sheet.  For 
each  element  in  the  sheet  we  first  insert  one  face  parallel  to 
the  boundary  face  and  all  the  edges  connecting  them  are  as¬ 
signed  with  zero  knot  intervals.  In  this  way,  only  the  bound¬ 
ary  nodes  have  non-zero  basis  function  values  on  the  solid 
T-spline  boundary  and  the  boundary  surface  will  be  defined 
by  the  boundary  nodes  only.  In  addition,  since  for  the  second 
boundary  layer  we  only  duplicate  the  nodes  originated  from 
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(a)  (b) 


(c) 


Fig.  4  The  boundary  layer  construction  operation  for  a  2D  example,  (a)  The  subdivision  result;  (b)  the  first  boundary  layer  for  the  subdivision 
boundary;  and  (c)  the  second  boundary  layer  for  the  input  boundary  surface. 


Fig.  5  The  boundary  layer  construction  result  for  the  Mouse  model,  (a)  The  T-mesh  after  boundary  layer  construction;  and  (b)  the  interior  of  the 
T-mesh.  The  pink,  blue  and  yellow  regions  are  the  subdivision  result,  the  first  boundary  layer  and  the  secondly  boundary  layer,  respectively. 


the  input  surface,  the  solid  T-mesh  boundary  is  exactly  the 
same  as  the  input  control  polygon.  In  the  following  steps  we 
will  not  modify  the  boundary  of  the  T-mesh.  Therefore,  we 
can  conclude  that  the  obtained  solid  T-spline  is  conformal 
to  the  input  T-spline  surface,  or  that  is,  the  input  boundary 
representation  is  preserved.  □ 


5  Irregular  Nodes  and  Solid  T-spline  Construction 

To  obtain  a  gap-free  T-mesh,  we  apply  templates  given  in 
[12]  to  the  interior  irregular  nodes.  There  are  two  kinds  of  ir¬ 
regular  nodes:  extraordinary  nodes  and  partial  extraordinary 
nodes.  Figure  6(a)  shows  the  template  for  a  partial  extraordi¬ 
nary  node,  in  which  the  magenta  edge  has  a  reflection  edge. 
Figure  6(b)  is  a  template  for  an  extraordinary  node.  We  in¬ 
sert  nodes  at  all  the  intersection  positions  for  each  element 
containing  irregular  nodes. 

Note  that,  in  the  obtained  solid  T-mesh,  the  layer  right 
below  the  boundary  only  contains  partial  extraordinary  nodes, 
and  the  edges  which  have  reflection  edges  with  respect  to 
these  partial  extraordinary  nodes  are  all  topologically  per¬ 
pendicular  to  the  T-mesh  boundary.  All  the  interior  extraor¬ 
dinary  nodes  are  at  least  two  elements  away  from  the  bound¬ 
ary  due  to  the  two  boundary  layers.  Hence,  we  will  not  in¬ 
sert  new  nodes  on  the  boundary  to  handle  extraordinary  and 


partial  extraordinary  nodes,  and  we  preserve  the  conformal 
boundary  property  in  this  step. 
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Fig.  6  Templates  for  a  partial  extraordinary  node  (a)  and  an  extraordi¬ 
nary  node  (b).  The  magenta  node  is  a  partial  extraordinary  node,  and 
the  red  node  is  an  extraordinary  node.  The  magenta  edge  is  the  edge 
with  a  reflection  edge  and  the  red  edges  have  zero  knot  interval. 

After  dealing  with  the  irregular  nodes,  the  local  knot 
vectors  for  each  node  are  inferred  and  for  each  domain  we 
detect  all  the  nodes  with  non-zero  basis  functions  and  use 
them  to  build  the  solid  T-spline  element.  A  rational  T-spline 
was  defined  in  [12]  whose  basis  functions  satisfy  partition  of 
unity  by  definition.  The  entire  solid  T-spline  model  is  built 
by  looping  over  all  the  local  domains.  For  the  constructed 
solid  T-spline,  the  continuity  on  the  boundary  is  the  same  as 
the  input,  and  the  continuity  is  C°  around  the  extraordinary 
nodes  or  partial  extraordinary  nodes,  and  C2  everywhere  else 
in  the  interior  region.  For  the  obtained  solid  T-spline,  we  ap¬ 
ply  Bezier  extraction  which  provides  a  finite  element  rep- 
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Table  1.  Statistics  of  all  the  tested  models. 


Model 

Input  surface 
nodes 

T-mesh 

nodes 

Interior 

extraordinary  nodes 

Interior  partial 
extraordinary  nodes 

Bezier 

elements 

Bezier  Jacobian 
(worst,  best) 

Time 

(s) 

Mouse 

(320,  164) 

1,429 

8 

44 

454 

(0.09,  1.00) 

1.9 

Boat 

(764,  420) 

2,852 

8 

56 

1,375 

(0.07,  1.00) 

5.2 

Propeller  Blade  (NURBS) 

(442,  256) 

1,763 

8 

44 

672 

(0.02,  1.00) 

3.2 

Propeller  Blade  (T-spline) 

(462,  272) 

2,083 

8 

64 

1,072 

(0.03,  1.00) 

4.8 

Propeller  Shaft 

(1,002,512) 

3,875 

8 

84 

1,556 

(0.05,  1.00) 

10.3 

Bunny 

(1,358,  1,083) 

7,537 

8 

85 

7,901 

(0.01,  1.00) 

32.1 

Fig.  7  Mouse  model,  (a)  The  constructed  solid  T-spline  and  T-mesh;  (c)  the  extracted  solid  Bezier  elements  with  some  elements  removed  to  show 
the  interior  mesh;  and  (c)  the  solid  T-spline. 


resentation  of  T-splines  for  isogeometric  analysis  [3,  8].  To 
construct  analysis- suitable  T-splines,  all  the  extracted  Bezier 
elements  must  have  a  positive  Jacobian  [14]. 


6  Results  and  Discussion 

We  have  applied  our  algorithm  to  several  models  (Figures  7- 
10).  Statistics  for  all  the  tested  models  are  listed  in  Table  1. 
The  results  were  computed  on  a  PC  equipped  with  an  Intel 
X3470  processor  and  8GB  main  memory. 

The  input  Mouse  and  Boat  models  were  designed  using 
the  commercial  CAD  software  Rhinoceros  and  the  inputs  are 
T-spline  surfaces.  Note  that  the  constructed  solid  T-spline 
can  preserve  all  the  sharp  features  in  the  input  model.  The 
Propeller  model  is  a  NURBS/T-spline  assembly  that  con¬ 
tains  three  parts:  two  blades  and  one  shaft.  The  input  model 
is  watertight  around  the  shared  boundary  of  the  three  parts, 
and  our  constructed  solid  T-spline  preserves  this  property 
exactly.  The  conformal  solid  T-spline  construction  can  be 
used  for  assembly  models  without  changing  its  mechanical 
fit  between  neighboring  parts.  The  input  of  the  bunny  model 
was  generated  using  earlier  work  in  [14]. 

For  all  these  models  the  obtained  solid  T-spline  has  ex¬ 
actly  the  same  spline  representation  as  the  input  surface, 
and  the  obtained  trivariate  T-spline  has  all  positive  Bezier 
Jacobians.  The  constructed  solid  T-spline  is  tricubic  and  C2- 
continuous  except  in  the  vicinity  of  partial  extraordinary  and 
extraordinary  nodes  in  the  interior.  On  the  boundary,  the 
solid  T-spline  has  the  same  continuity  as  the  input  NURBS/T- 
spline  surface. 


7  Conclusions 

We  presented  a  novel  method  to  construct  conformal  solid 
T-splines  for  genus-zero  geometry  from  a  boundary  T-spline 
surface.  Based  on  the  boundary  layer  construction  technique, 
the  solid  T-spline  has  the  same  spline  representation  as  the 
input  spline  surface  and  can  preserve  sharp  features  in  the  in¬ 
put  model.  Our  method  is  efficient  and  robust.  In  the  future, 
we  intend  to  extend  the  algorithm  to  more  general  geome¬ 
tries  with  arbitrary  topology  and  construct  analysis- suitable 
solid  T-splines. 
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Fig.  9  Propeller  model,  (a)  The  input  T- spline  boundary  surface  with  the  control  mesh;  (b)  the  constructed  solid  T- spline  and  T-mesh;  (c)  the 
extracted  solid  Bezier  elements  with  some  elements  removed  to  show  the  interior  mesh;  (d)  the  extracted  solid  Bezier  elements;  and  (e)  the 
constructed  solid  T-spline. 
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Fig.  10  Bunny  model,  (a)  The  input  T-spline  boundary  surface  with  the  control  mesh;  (b)  the  constructed  solid  T-spline  and  T-mesh;  (c)  the 
extracted  solid  Bezier  elements;  (d)  the  extracted  solid  Bezier  elements  with  some  elements  removed  to  show  the  interior  mesh;  and  (e)  the 
constructed  solid  T-spline. 
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